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Abstract. Effects of randomness on the spin-i and 1 antiferromagnetic Heisen- 
berg chains are studied using the quantum Monte Carlo method with the 
continuous-time loop algorithm. We precisely calculated the uniform suscep- 
tibility string order parameter, spatial and temporal correlation length, and 
the dynamical exponent, and obtained a phase diagram. The generalization of 
the continuous-time loop algorithm for the systems with higher-5 spins is also 
presented. 

1. Introduction 

Low-energy properties of one-dimensional quantum Heisenberg antiferromag- 
nets have been of great interest for many years. Particularly, effects of random- 
ness on ground states with strong quantum fluctuations, such as the so-called 
Haldane state, have been studied extensively in recent years [1-5]. 

In the present work, we investigated S = i and S = 1 random antiferro- 
magnetic Heisenberg chains by means of the quantum Monte Carlo method 
with the continuous-time loop algorithm [6-9]. To this end, we first generalize 
this algorithm to systems with S — 1 or higher. This quite powerful simulation 
technique enables us precise calculations of thermodynamic quantities, such as 
the uniform susceptibility, etc. 

2. Model and Numerical Method 

We consider the spin-5 1 antiferromagnetic Heisenberg chain of L spins with 
random nearest-neighbor interaction: 

L 

H = Y,Ji{Si -Si+i), (1) 

i=l 

where Si is a spin operator at site i and Sf = S(S + 1). A periodic boundary 
condition is applied. All of Jj's are taken to be positive (antiferromagnetic). 

It is well known that traditional world-line Monte Carlo methods [6] suffer 
from long auto-correlation time at low temperatures. This difficulty is solved by 
the loop algorithm [7,8]. In addition, one can completely eliminate the system- 
atic error beforehand, which is due to the discreteness in the imaginary-time 
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direction, by adopting the recently proposed continuous-time algorithm [9]. 
The continuous-time loop algorithm, which was originally proposed for a spin- 

1 system, can be generalized in a straightforward way to the spin-1 or higher 
systems as shown below [10]. 

First, we introduce spin-^ representations of the partition function of the 
Hamiltonian (|l|): 

Z = Y, <«l exp[-/3H]P |n) , (2) 
M 

where 

L 2S j 

n = ^2^2 ~j( a w ' ^i+i.f) • ( 3 ) 

i— 1 1 

Each Si is now represented as a sum of 25* Pauli operators (c^i, <Ti t 2, 
Ci,2s) and the trace in Eq. @ is taken over the complete basis {n} of the 
space spanned by the 2SL Pauli operators. The projection operator P = Y[ i Pi 
projects out unphysical states, which do not appear in the original phase space 
of dimension (2S+1) L . After performing the Suzuki- Trotter decomposition, we 
obtained the partition function of the (l+l)-dimcnsional classical Ising model 
with four-body weights of exp[— /3,/j(7 ijM <7j +lj „/4m] and additional (4S')-body 
weights of Pi at the boundaries in the imaginary-time direction. 

Note that each four-body weight is completely equivalent to that appears in 
the world- line representation of the spin-^ system. Therefore, in constructing 
a loop algorithm, assignment probabilities of graphs for these weights remain 
unchanged. For a boundary (4S')-body weight, we consider (25)! types of 
graphs, each of which consists of 2S edges connecting one of | spins at r = /3 
with a j spin at r = one by one. A graph is called 'compatible,' if every 
pair of spins connected by an edge have a same direction. A graph is selected 
out of such compatible graphs with equal probability. After the assignment 
procedure, loops are identified to be nipped independently with probability \. 

The continuous-time limit (so-called Trotter limit) of the above procedure 
clearly exists. This is the most significant advantage of the present algorithm 
when one makes a comparison with the discrete version of the general-5 1 loop 
algorithm, proposed by Kawashima and Gubernatis [12]. The present method 
also works with the XXZ model, systems in higher dimensions, and even sys- 
tems with mixed magnitude of spins without any changes on the algorithm. 

We have calculated the spin-spin correlation function of the clean antifer- 
romagnetic Heisenberg chain with S — |, 1, |. The correlation length in 
the zero temperature limit is estimated as £ ~ 6.02 and 49.6 for S = 1 and 

2 cases, respectively, while we observed power-law divergence with the inverse 
temperature in the cases with half-integer spins. The estimates for S = 1 and 
2 agree quite well with the recent results by the density matrix renormalization 
group method [13,14]. 
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Figure 1: Uniform susceptibility in the RS phase of the S = \ random antifer- 
romagnetic Heisenberg chain with bond distribution (W, m)=(1,0) and (2,0). 



3. Results 

We consider the uniformly-distributed bond randomness with dimerization: 
Ji = 1 + (-!)*« + Si (4) 



with 



P(*) = j 1/W (5) 
otherwise . 



The parameters W and u control the strength of the randomness and the 
dimerization, respectively. For |u| + 5- < 1, all Jj's are antiferromagnetic. 

3.1. S = i random Heisenberg chain 

The clean (i.e. W = and u = 0) S = | antiferromagnetic Heisenberg chain 
has no excitation gap. By introduction of infinitesimal randomness, the system 
is driven to the random singlet (RS) phase. In this phase, it is shown that 
there is also no excitation gap, but the antiferromagnetic correlation of spins 
decays with a different exponent from that in the clean case [1,2]. 

In Fig. |l|, we show the uniform susceptibility in the cases with (W, u) = (1, 0) 
and (2,0). The random average was taken over 3000 samples for the lowest 
temperature (T = 0.001) in the calculation. The uniform susceptibility diverges 
as the temperature goes down, which makes a sharp contrast with the clean 
(W — 0) case. As seen in Fig. [ij (x^ 1 ) -1 ^ 2 plotted versus logT lies on a straight 
line at low temperatures in both cases. This implies the uniform susceptibility 
behaves as \ ~ l/T(logT) 2 , which is consistent with the prediction of the 
real-space renormalization group (RSRG) theory [2]. 
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Figure 2: Exponent 7, spin gap A, and square root of the inverse spatial 
correlation length £ -1 / 2 of the S — i dimerized random antiferromagnetic 
Heisenberg chain with u + -5- = 1 . 



3.2. S = i random chain with dimerization 

If a weak dimerization u is enforced with u+ fixed to 1, the spin correlation 
becomes short ranged immediately, while the spin gap remains closed until 
u ~ 0.27 (Fig. ||). This phase is referred to as the random dimer (RD) phase, 
and a realization of the quantum Griffiths phase [2,3]. In the RD phase, the 
exponent 7 of the uniform susceptibility (x ~ T^ 1 ) varies continuously and 
diverges to the negative infinity at the end of the RD phase. 

The string order parameter O str exists in the dimer phase and also in the 
RD phase. We found near the RS point (u — 0), O s tr is well scaled by a 
logarithmic scaling form: 

O str ~ L-W/'fiuL 1 /", L x l v I logT) (6) 

with v = 2 and f3 — 0.382 as predicted by the RSRG arguments [2,3]. This 
implies the dynamical exponent z is actually infinite in the RS phase. 

3.3. S = 1 random Heisenberg chain 

Finally, we consider the S = 1 chain without dimerization. We found almost 
similar behavior of the uniform susceptibility as in the case of the S = i 
dimerized case discussed above. For small randomness W, the system has 
the so-called Haldane gap. As the strength of the randomness increases, the 
exponential behavior of the susceptibility changes to the power-law behavior 
at W - 1.5. 

For the extremely wide bond distribution (W = 2), we have confirmed that 
the uniform susceptibility obeys the same asymptotic form \ ~ l/T(\ogT) 2 as 
in the RS phase of the S = % case. 
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4. Summary and Discussions 

Results of the precise quantum Monte Carlo simulation of the random anti- 
ferromagnetic Heisenberg chain have been presented. We have obtained the 
phase boundary between the RD and the dimer phases on« + y = 1- We 
expect there would be a phase boundary line connecting (W, u) — (0, 0) and 
(1.46,0.27) in the W~u plane. 

For | it | < the bond distributions of odd and even bonds overlap with 
each other. Therefore, it is naturally expected that one could find a very long 
and almost clean segment with arbitrary small excitation gap. However, the 
present result clearly shows that the spin gap can survive even in such regions. 
Its proper interpretation remains as an open question. 

For the 5 = 1 chain, similar behavior of the uniform susceptibility as in 
the case of S = \ dimerized chain has been found. Detailed analysis (critical 
strength of randomness, etc.) for this system is now being proceeded and will 
be presented elsewhere. 
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